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Motivated by the experimental observations of resonant tunnelings in the systems with 
half-integer spin, such as V15 and Mn4, we study the mechanism of adiabatic change of 
the magnetization in systems with the time-reversal symmetry. Within the time-reversal 
symmetric models, effects of several types of perturbations are investigated. Although tun- 
neling between the ground states is suppressed in a simple Kramers doublet, we show that 
the nonadiabatic transition governed by the Landau-Zener-Stiickelberg mechanism occurs 
in many cases due to the additional degeneracy of the ground state. We also found more 
general cases where LZS mechanism can not be applied directly even the system shows a 
kind of adiabatic change of the magnetization. 

§1. Introduction 

Recently the dynamics of nanoscale molecular magnets has attracted much inter- 
est. In particular, the magnetization processes in some of molecules show step-wise 
structures of hysteresis in the sweeping magnetic field. This phenomenon is at- 
tributed to the quantum mechanical resonance at level crossing points in discrete 
energy structure of the small magnets, and is called resonant tunneling. The low- 
energy magnetic properties of such molecules is usually given by a model of a single 
large spin. For example, the system of Mni2 ^ or Fe§ 2 **' 3 ** is represented by 

H = -D{S Z ) 2 - HS Z , S = 10, D>0. (1-1) 

Level-mixing interactions such as 

Hi = C((S+) 4 + (S-)*) or Hi = E((Sf ) 2 - (Sf) 2 ). (1-2) 

cause hybridization of the degenerate states of different values of the magnetization, 
and the resonance between the two states takes place, i.e., the resonance tunneling. 
The hybridization causes an energy gap, i.e., the tunneling gap. In particular, in the 
case of integer spin S, the hybridization occurs between two-fold degenerate ground 
states with the total magnetization M = ±S at H = 0. This tunneling phenomena 
has been studied by a kind of path-integral method 4 ) which has provided various 
important informations. There the tunneling probability is expressed by a form of 
path integral representation, e.g., 

P= (vr|e-^|0), (1-3) 
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where H = K X (S X ) 2 + K y S y , K x > K y > and \tt) and |0) denote the states 
directing ±z directions (the easy axis), respectively. 

However, in this picture of a single large spin model, the tunneling in half-integer 
spin particle is known to be suppressed due to the interference of the Berry phase 
5 ) " 7 ) . This fact is originated from the time-reversal symmetry and it corresponds to 
the Kramers theorem that all the energy levels must be at least doubly degenerate 
in half-integer spin systems. In the context of eq.(T3), the contributions from the 
pair of the two paths related by the time-reversal operation cancel with each other 
and 

(7r|e^ w |0) =0, (1-4) 

when S = half-integer, from which one may conclude the absence of the tunneling 
gap. 

Generally, systems have the time-reversal symmetry when the hamiltonian of 
the system consists of products of even number of spins, e.g., the usual two spin 
interactions, etc. Therefore, in systems of odd number of spins, the eigenstates must 
be at least doubly degenerate at H = 0, i.e. the Kramers theorem. However, in many 
half-integer spin systems, e.g. Vis(5 = 3/2) 8 )' 9 ) and Mn^S 1 = 9/2) 10 ), temperature 
independent relaxation phenomena have been found, which suggests the existence of 
a tunneling splitting in half-integer systems. 

In this paper, we consider this tunneling problem in time-reversal symmetric 
systems. The key point is the fact that the realistic molecule consists of many 
atoms and many degrees of freedom exist. There the magnetic property can not be 
simply expressed by a single large spin. We find that the system can show tunneling 
behavior making use of additional degrees of freedom, even if each level are the 
Kramers doublet. 

The existence of the tunneling gap causes an avoided level crossing structure of 
the energy level as a function of the magnetic field H 

Hz = -HJ2S?, (1-5) 

i 

and we expect the adiabatic motion of the magnetization when we sweep the mag- 
netic field. As time-dependent phenomena, the magnetization adiabatically follows 
the ground state value if the field-sweeping is slow. There the sign of the magneti- 
zation changes near H = 0, which is the adiabatic motion. If we sweep H fast, then 
generally the nonadiabatic transition occurs 11 )" 21 ), and the transition probability of 
the nonadiabatic transition was given by Landau, Zener and Stiickelberg n ) " 13 ) as 

p = exp(-^), (1-6) 

where v is the sweeping rate of the field dH/dt, and AE is the energy gap at the 
avoided level crossing structure. This transition rate plays an important role in the 
resonant tunneling phenomena in nanoscale molecular magnets such as Mni2 and 
F eg 3),i7)-2i)^ j n ^ e present paper we point out that the various types of avoided 
level crossing structures are possible in half-integer spin systems which are time- 
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reversal symmetric at H = 0, and also we study properties of adiabatic changes of 
the magnetization. 

§2. Effects of magnetization-nonconserving matrix elements 

In this paper, we mainly consider the Heisenberg model with odd number of 
spins as the unperturbed system. In order to change the value of the magnetization, 
the system must contain some terms which violates the conservation of the magne- 
tization. That is, the hamiltonian must not commute with the total magnetization. 
We will see that some of such systems have the avoided level crossing structures. 

Here we study an antiferromagnetic Heisenberg chain for Ho 

(ij) 

with various types of perturbations Hi of the form 

H. »tV s J- (2-2) 

ij 

First we study a simple case of three spins (N = 3), and more complicated cases will 
be studied later. The states of the pure Heisenberg model (ccjj = 0) are classified 
by the total spin S, that is, the four states of S = 3/2 and two sets of two states 
of S = 1/2 for the present case of N = 3. In Fig. 1(a), we show the energy level 
structure of the system for a = 0.0 as a reference. Here we see six lines but the levels 
crossing at low energy consist of two degenerate states. That is, there are four states 
of S = 1/2. This degeneracy for the S = 1/2 states appears due to the rotational 
symmetry. 

When we apply the perturbation (2-2), the total spin is not a good quantum 
number. However, it gives still good description of the states as far as the perturba- 
tion is small. Thus hereafter we also use S to denote the states. 

Here we study the following five categories of the perturbation: 



Model 


Rotational symmetry 


Reflection symmetry 


I 


O «12 = Q23 ' ' ' = ONI = a 


O aij = aji 


II 


X 


O Qi2 = Q21 = a, others = 


III 


o 


x aij = —aji = a 


IV 


X 


x ai2 = — Q21 = a, others = 


V 


X 


x ai2 7^ Q2i, others = 



Table 1 Types of perturbations 



2.1. Model I 

In Fig. 1(b) we show the energy level structure for Model I with a = 0.2. There 
we find that the hybridization does not occur at lower crossing point while it occurs 
at higher crossing point. In this model the perturbation is symmetric for both the 
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Fig. 1. (a) Energy structure of the pure Heisenberg model of N 
a = 0.2 



3, and (b) of the Model I with 



(2-3) 



rotation and reflection, and does not affect the symmetry of the wavefunctions. Here 
we have 

Wi|3/2>= |3/2) 
Wi|l/2>= |l/2), 

where 1 3/2) , and |l/2) denote the Hilbert space of the wave functions with S = 3/2 
and 1/2, respectively. In particular, in the space of S = 1/2 the symmetric operation 
H\ does not have matrix elements 



Wi|l/2) = 0. 



(2-4) 



Therefore the degenerate crossing structure of S = 1/2 states in Fig. 1(a) is preserved. 
On the other hand, in the space of S = 3/2, Hi has matrix elements between 
M = 3/2 and 1/2, and between M = —3/2 and —1/2, where M is the z component 
of the magnetization. The absence of matrix elements between M = 1/2 and —1/2 
is due to an peculiar interference of operations of this type of Hi (see Appendix). 
Thus, avoided level crossing structures are formed only between the states M = 3/2 
and 1/2 and between M = -3/2 and -1/2. 

Thus if we consider the antiferromagnetic model of the type of Model I, the 
adiabatic transition can not occur between the M = 1/2 and —1/2 states as we 
expected from the single spin model. 

In the present paper we will not mention for the ferromagnetic model in detail. 
However, it should be noted that the property of the ferromagnetic case is simply 
obtained by reversing the energy axis. In the present model, adiabatic transitions 
occur between M = -3/2 and -1/2, and M = 1/2 and 3/2. 



2.2. Model II 



Model II of a = 0.2 shows a similar structure to that in Fig. 1(b) as shown 
in Fig. 2(a). However due to the lack of the rotational symmetry, it shows more 
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Fig. 2. (a) Energy structure of the Model II of N = 3 with a = 0.2 (b) with a = 0.8 

complicated structure as shown in Fig. 2(b), where we use a = 0.8 to emphasize 
the crossing structure. This model has the reflection symmetry with respect to the 
exchange IZu of the sites 1 and 2. Thus the states are classified as 

where |+) and |— ) denote the symmetric and antisymmetric states with respect to 
the exchange IZu- The four states in S = 1/2 are separated into sets of |+) and |— ). 
The antisymmetric states |— ) do not contain all up state or all down state, and are 
kept in strict S = 1/2 space. Thus they are not affected by Hi and degenerate at 
H = 0. These two states of +) in S = 1/2 contain a little components of the all up 
and down state and no more are in S = 1/2 space strictly. However the two states 
of |+) also must degenerate at H = because of the time-reversal symmetry. In this 
model, the adiabatic transition can not occur between the M = 1/2 and —1/2 states 
when the field crosses zero similarly to the previous model. The state of S = 3/2 
also forms complicated avoided structure. 

In Fig. 3(a), we show the magnetization of the states of four low energies as 
functions of the field for the case of a = 0.8. The bold solid line, bold dotted-line, 
dashed-line and thin solid line denote the magnetization of the levels 1,2,3 and 4, 
respectively. When the levels cross the order of the lines changes. Therefore the types 
of lines change discontinuously. However, the adiabatic motion of magnetization 
is given by smooth continuation of a level. The adiabatic magnetization starting 
with the ground state at H = — 5 is shown by big circles. It shows a strange field 
dependence. That is, around H ~ 2 the magnetization has a peak, where the 
coefficient of the all up state becomes large. 

When we sweep the field, 



H(t) = -H + vt, 



(2-6) 
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the evolution of the state is expressed by 

\V(t)) = exp t (-ih H(s)ds^J \V(t)), H(t) = n-H(t)J2Sf, (2-7) 

where exp t means the time ordered exponential. In Fig. 3(b), the magnetization 
process with sweeping field with velocity v = 0.02, 0.2 and 2.0, are shown. In the 
case of v = 0.02 the magnetization in the time dependent process reproduces well 
the adiabatic one. 




Fig. 3. (a) Adiabatic magnetization process of the Model II with a = 0.8, (b) Magnetization process 
with sweeping field v — 0.02, 0.2 and 2.0, shown by the solid line, dotted-line and thin solid line, 
respectively. 

Here, it should be noted that the magnetization process shows a non-monotonic 
dependence on the field, which violates the thermodynamical stability. That is, the 
Zeeman energy (1-5) decreases when the field increases. Namely, if we increase the 
magnetic field the magnetic field receives some energy but not gives the energy as in 
the normal cases. This thermodynamically unstable behavior comes from the fact 
that this adiabatic state is no more the ground state. In the ground state the Zeeman 
energy must increase monotonically with H. In the present case, the ground state 
for H < adiabatically becomes an excited state after the crossing point due to the 
crossing. It would be interesting to find phenomenon corresponding to the present 
observation and to make use this peculiar behavior of the magnetization. 

2.3. Model III 

In Fig. 4(a), we show the energy level structure for Model III with a = 0.2. 
Using a similar argument for Model I, 

Wi|3/2)=0, (2-8) 

and thus we find that the structure of S = 3/2 states is preserved. On the other hand, 
Hi can change the magnetization in the space of S = 1/2. Here the states of S = 1/2 
with different magnetizations M = ±1/2 are hybridized and form two sets of avoided 
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Fig. 4. (a) Energy structure and (b) ground state magnetization, of the Model III of N = 3 with 
a = 0.2 



level crossing structure. The magnetization process of the ground state is shown in 
Fig. 4(b). They are completely degenerate due to the rotational symmetry. The 
adiabatic transition between the different sets of avoided cross structure is prohibited. 
Let us take a ground state \a) at H = — Hq in one set of avoided cross structure 
as an initial state. We have another ground state \b) in the other set which is 
orthogonal to \a). The orthogonality does not change by the sweeping operation 
because of the symmetry. For example if we begin with a state whose eigenvalue is 
g 4j7r/3 £ or Qne third rotation, the state with the eigenvalue e 8l7r / 3 is not mixed at all. 
Thus we can consider the time evolution for both states separately. In each branch, 
the nonadiabatic transition probability is exactly given by the equation (T6), and 
thus the total change of the magnetization for any combination is given by the LZS 
formula. 

2.4. Model IV 

In Fig. 5(a), we show the energy level structure for Model IV with a = 0.2. 
There we find a similar structure to that of Model III. However, because of the lack 
of the symmetry, the four levels of low energies are no more degenerate and if we 
look the state carefully they have different field dependences (see the inset). In order 
to emphasize the structure, the energy structure of a = 0.8 is also shown in Fig. 
5(b). Although this model has no geometrical symmetry, the model has a peculiar 
symmetry due to the special interference in the triangle lattice. That is, the all up 

state | +++) can not reach to the all down state | ). Thus the states is classified 

into the following two categories 

|A)= a| + ++) + 6|+) + c|-) e 

|B)= a'\ ) + 6'|-)o + c'|+) e , l/ ' yj 

where |+) e is a symmetric state of M = 1/2 with respect to the exchange of the 
sites 1 and 2, and |— ) Q is the antisymmetric state of M = —1/2, and so on. In 
the inset of 5(a), the open circles denote the level of |A) and the closed circles |B). 
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Fig. 5. (a) Energy structure of the Model IV with a = 0.2, (b) with a = 0.8 



Avoided crossing structure is formed in each category. Thus there are no mixing 
between different sets of the avoided crossing structure. Because of the Kramers 
theorem, they cross at H = and recover the degeneracy. For any combination 
of the initial populations of the two sets, the nonadiabatic transition probability in 
each set is given by (1-6), and the change of the total magnetization follows the LZS 
mechanism. In this model, if we sweep the field in a wide range, e.g., from —4 to 
4, the magnetization adiabatically changes from M = — 3/2 to M = 1/2, but not 
to 3/2. This asymmetry of the adiabatic change is one of the characteristics of the 
present model. 

Up to the model IV, the system has some symmetry with respect to the exchange 
of spin 1 and 2. Now we try to remove this symmetry. Now we put only aSfS^, but 
there we found essentially the same structure to that of Model IV. We found that 
there is still some symmetry which causes the decoupling of the space into 2+2. Even 
if we add a term fiS^S^ in order to reduce possible symmetry, we find that there still 
exists two sets of states. These separations are due to rather special interference on 
the triangle structure. 

2.5. Model V 

Finally we put a\2 = 0.2 and a 2 \ = —1.0. In Fig. 6, we show the energy level 
structure for this case, which is similar to the previous one. But in this model we 
found that there is no separation into 2 + 2 subspaces any more. The overlaps of the 
ground state at H = —0.3 (\G(H = —0.3))) and the eigenstates at H {\4>k{H)) at 
H) of four lowest eigenvalues are measured by the quantity 



x(k) = \(G(H = -0.3)\MH))\ 2 , k = 4, (2-10) 
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as shown in Figs. 7. In Fig. 7(a), we 
show the overlaps in a case of ct\2 = 0.2 
and «2i = where the separation into 
2+2 occurs, and in Fig. 7(b) the present 
case (a\2 = 0.2 and ol<i\ = —1.0). Here 
the data for k = l,---,4 are shown 
by bold-solid, bold-dotted, dashed and 
thin-solid lines, respectively. While the 
set of states of k = 1 and 3 and that 
of k = 2 and 4 do not mix at all in 
Fig. 7(a), all the four states have some 
overlap with the state \G(H = —0.3)} 
in Fig. 7(b). The bold line in Fig. 7(b) 
shows very small but nonzero values at 
H > 0. The wave functions consist of 
all the states 

\V(t)) =a\ + ++) + 6|+) + c|-) e + o'| ) + b'\-) Q + c'|+) e . (2-11) 

Thus, the simple LZS transition mechanism does not work in this model. Dynamical 
property of this model will be studied in the next section. 




/ 

/ 
/ 

/ 



-1 o 1 

H 



Fig. 7. (a) Overlap x(k) for an = 0.2 and 021 = and (b) for 012 = 0.2 and ai\ = —1.0 
2.6. Summary 

Here we study the simple case of N = 3, but we can learn general properties for 
higher spins. That is, although the time-reversal symmetry prohibits the tunneling 
between the Kramers doublets at H = 0, it does not means the absence of the 
adiabatic change of magnetization when the field swept at H = 0. If the system 
has additional degree of freedom and more than two states are degenerate, e.g. 4- 
fold degeneracy, at H = in the unperturbed system Ho, the avoided level crossing 
structure is formed with respect to the field as in the cases of Model III, IV, and V. 
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If the system has additional symmetries, then the avoided level crossing structures 
are classified by the symmetry and they are independent each other, where the LZS 
transition mechanism works as in the two-level system. On the other hand, in the 
absence of additional symmetry, the degenerate states forms an irreducible space 
except H = 0. There the simple LZS transition mechanism does not work. In the 
next section we study the magnetization change under the sweeping field in the latter 
case. 

§3. Quasi-Landau-Zener-Stiickelberg Transition 




Fig. 8. (a) Magnetization process in the sweeping speed v = 0.04, 0.08, 0.12, and 0.4 from the 
top to the bottom. ai2 = 0.02, ai\ = —0.1 (b) The overlaps with the k-th adiabatic state 
Pk(H) for k=l,2,3 and 4 by the bold solid, bold dotted, dashed and solid lines, respectively. 
ai2 = 0.02,a2i = -0.1, v = 0.04. 

Now we demonstrate the time dependent process with sweeping field where H 
is swept from — Hq to Hq. 

|!^(t f )) = exp t [i £ ! H(s)ds S j \G(H = -H )), (3-1) 

where \G(H)) is the ground state for H, and H(t[) = —Hq and H{t{) = Hq. In 
Fig. 8(a) we show the magnetization for the case of au = 0.02 and «2i = —0.1 
(Model V). Here we take small values of a's in order to facilitate the observation of 
the nonadiabatic transition. The speed v of the field-sweep is 0.02, 0.04, 0.08, and 
0.2 from the top to the bottom, respectively. As we saw in the previous section, the 
simple LZS mechanism can not be applied for this model because all the four low 
energy states contribute to this transition. However we find similar ^-dependence 
of the magnetization change to that in the LZS mechanism. In order to study the 
detail of the transition we investigate the population of the fc-th eigenstate at the 
field H by 

Pk (H) = \(*(t)\MH))\ 2 , H = H(t). (3-2) 
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In Fig. 8(b), the change of the populations are shown, where we find the ground 
state at large negative field is scattered to all the four state and they change with 
the field. This is a characteristic feature of the present model in contrast with the 
other cases: e.g., in the case of «2i = 0, where the four states are separated into two 
sets and the ground state at high negative field is scattered only to the second and 
forth states. 



1 



Po( v ) 

0.5 




0.2 0.4 
V 

Fig. 9. Nonadiabatic transition probabilities in the Model V with ai2 = — a2i/5 = 0.03, 0.02, and 
0.01, from the top to the bottom. The solid curves represents the value estimated by the LZS 
formula (1-6). 

However, if we sum up the populations of the first two states, 

po = mtf)\MH))\ 2 + mt f )\MH))\ 2 , (3-3) 

the sum almost satisfies the dependence of LZS transition probability (1-6) as shown 
in Fig. 9 where 

Po = l-p. (3-4) 

Here we use the energy gap at H = as AE. The data for the system of au = 
— Q2i/5 = 0.03, 0.02, and 0.01 are shown by circle, square, and triangular, respec- 
tively. The lines denote the corresponding dependences due to the relation (1-6). 

This is a problem of nonadiabatic transition of four levels. But the mechanism 
of the present model seems to belong to a different type from that of the bow-tie 
model 23 ) . The present transition may be attributed to some multiple crossings 24 )' 25 ) , 
and then have some relation with the Brundobler and Elser hypothesis 26 ) . The 
mechanism of this quasi-LZS behavior of the sum of the population would be an 
interesting problem in the future. 

From a view point of the present study, the resonant tunneling from —1/2 to 1/2 
at H = in molecular magnets with half-inter spin such as V15 does not contradict 
with general theory of the symmetry. Here it would be interesting to point out the 
following fact. If the field is swept slowly from a large negative value where the mag- 
netization is —3/2, the magnetization change adiabatically as —3/2 — > —1/2 — > 1/2, 
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but not to 3/2 as shown in Fig. 10(a), although the ground state magnetization 
process is, of course, symmetric (Fig. 10(b)). We found this behavior exactly in the 
Model IV. However, in the present model the wave functions contain both compo- 
nents of | + ++} and | ). Therefore this asymmetry of the adiabatic magne- 
tization seems to occur widely in model on the triangle. If the system essentially 




Fig. 10. (a) Adiabatic magnetization process (v — 0.002) and (b) Magnetization process of the 
ground state, of Model V for ai2 = 0.02 and cti\ = —0.1 

consists of three spins, this asymmetric adiabatic process would be observed in a 
slow sweeping field at low temperatures. However, if other spins strongly contribute 
to the interaction, the states of 3/2 could be also degenerate. In such cases also the 
transition from 1/2 to 3/2 would be allowed. It would be an interesting problem 
to study the magnetization change of V15 in a slowly sweeping field , i.e., adiabatic 
change of the field. This enables us to determine whether only three spins in the 
middle of the molecule mainly contribute to the interaction, or the whole 15 spins 
contribute. We hope that some experimental study on this problem will be done, 
although the thermal effects would easily smear out the pure quantum selection rule. 

§4. Higher Spin Cases 

So far we studied magnetic proper- 
ties the crossings which are mainly con- 
sists of S = 1/2 states. Here let us 
study the cases of higher spins. As an 
example, we study a system of five spins 
shown in Fig. 11, where the bold line 
denotes J, and the thin line denotes J'. 
The perturbation (2-2) is set between 
the sites 1 and 2 only with a\2 = 01. 
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If we set J = J' = 1 and a = 0, the 
ground state at H = is given by 4-fold 

degenerate S = 1/2 states, and the energy structure is similar to Fig. 1(a). If we 
put a = 0.2, the energy structure is similar to Fig. 5(a), and we find two sets of 
independent avoided level crossing structure again. Thus, it is expected that the LZS 
transition works in a wide range of systems when the ground state is in S = 1/2. 

In order to have the ground state with a higher spin in antiferromagnetically 
interacting system, we set J = 1, J' = 0.2. For a = 0, the ground state at H = is 
given by 4-fold degenerate S = 3/2 states. The energy structures for a = and 0.2 
are given in Fig. 12(a) and (b), respectively. The zoom-up around H = is shown 
in Fig. 13(a), where we find a complicated structure. 




Fig. 12. (a) Energy structure of the pure Heisenberg model and (b) of the Model V with a = 0.2 

The overlaps among all the four states are nonzero. In Fig. 13(b) the overlap 
functions are shown as in the same manner as eq.(2-10), where the solid line, bold 
dotted-line, dotted-line, dashed-line, dot-dashed line, and dot-dot-dashed line denote 
x(k), k = 1, • • • , 6, respectively. Here we do not expect simple LZS transition. 

Although we will report elsewhere detail properties of this non-LZS type adia- 
batic transition of high spin states, we show the adiabatic magnetization processes 
of the present model in Fig. 14(a). If we sweep the field from a large negative value, 
the adiabatic magnetization is given by the bold solid line for H < 0. This bold solid 
line jumps at H = 0. But the ground state for H < adiabatically continues to the 
second level, and there the adiabatic magnetization is given by the bold dotted-line. 
The second level forms an avoided level crossing with the third level near H = 0.02. 
If the nonadiabatic transition occurs at this point, the magnetization is given by the 
dashed-line instead of the bold dotted-line after this point. Thus we expect that if 
v « 1 then the magnetization curve is given by the bold dotted-line, and if v » 1 
then by the dashed-line. If the field is swept with a finite speed, the magnetization 
is located between the bold dotted-line and the dashed-line depending on the speed. 
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In Fig. 14(b), the magnetization change with sweeping velocity v = 0.0002 is shown, 
which is indeed between the dotted-line and the dashed-line. 

Here we adopt a peculiar example for higher spin case. The property of the 
crossing would depend on the types of perturbation. It should be also noted here 
that so far we study only isotropic model for Ho. However, when we study high spin 
cases, the uniaxial anisotropy (1-1) plays an important role, too. Thus, we expect a 
variety of types of adiabatic transition of the magnetization in high spin cases, which 
will be reported elsewhere. 




Fig. 14. (a)Magnetization of low energy levels as a function of the field H, (b) Magnetization change 
with sweeping field v = 0.0002 
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§5. Summary 



We have studied the resonant tunneling phenomena in half-integer spin systems 
with the time-reversal symmetry. If the system has only single degree of freedom, 
the time-reversal symmetry prohibits the adiabatic change of the magnetization as 
has been pointed out. However, if the system consists of many spins and there exist 
another degree of freedom, various types of adiabatic change of the magnetization 
occurs. 

In the present paper, we found degenerate (or quasi-degenerate) avoided level 
crossing structures in Model III (or IV, V), where the Landau-Zener-Stiickelberg 
mechanics works. Appearance of avoided cross structures in half-integer spin sys- 
tems does not contradict with the Kramers theorem in all cases because the states 
belonging to different avoided level crossing structures degenerate at H = 0. 

The simple LZS mechanism of the nonadiabatic transition does not work in some 
cases where more than two levels are involved in a nonadiabatic transition (Model V). 
However an interesting sum rule seems to hold, namely when we add the transition 
probabilities to the two neighboring states, the LZS formula gives an accurate fit 
to the numerically obtained data. This issue deserves further studies as a possible 
extension of the LZS theory. 

We would like to point out that the adiabatic magnetization does not necessarily 
coincide with the magnetization of the ground state even the adiabatic change of the 
sign of the magnetization occurs at H = 0. When we sweep the field from a large 
negative value where the magnetization is —S, we could not find any case where the 
magnetization reaches to S when the field becomes a large positive value. Although 
the naive extension of (1-4) does not hold when we sweep the field, i.e., 



where the H(0) = —Hq < and H(t) = Hq, the present observation may suggest 
that the above matrix elements is zero in the limit Hq — > oo. This problem will be 
studied in more general cases. 

As we saw in Model I and II, if the system has some symmetry and there are 
only two state in the same symmetry, the levels simply cross at H = due to the 
time-reversal symmetry. Even in such case the adiabatic change of the state can 
provide nontrivial behavior at H / as we saw in Model II. The thermodynamical 
unstable state would cause interesting phenomena when the system is coupled with 
the dissipative environments. 
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Appendix A 

Special interference on the triangle lattice 

Let us consider Hi of Model I 

Hi = a(SfS z 2 + SfSi + Si Si + Sf Sf + Sf Sf + S z 3 Sf). (A-l) 
If we consider the total spin operators 

QX _ QX , QX , r>2/ _ Cf I QV I C2/ Q2 _ Q2 , QZ , QZ /A r>N 

°tot 1 ' 2 "l" 3 ' J tot 1 ' 2 3 ' J tot 1 "i - 2 ' 3' V A Z J 

Hi is expressed as 

^1 = ^ ('S'tot'S'tot + ^tot'S'tot) = a ^tot'S'tot + 2^*0*^) ' 

If we apply Hi to the all-up state | + ++), we have 

Hi\ + ++) = a (Sfa x § + iS^)) I + ++} 

= « (f^tot - R ot ) I + ++) (A-4) 
= f 'S'totl + ++)• 

This is a symmetric state in M = 1/2. If we apply Hi further, we have 

H\\ + ++) =i (SZ* x \ + iS^)) S t - ot | + ++) 

= ^ (^fot + £Sg*) S^tl + ++) 
= ^+5^1 + ++)oc| + ++). 

Thus, the space of M = —1/2 can not be reached from | + ++) by Hi- 

The present argument is directly extended to general N spin system with inter- 
actions between all spin pairs. There we find that Hi in (A-3) has no matrix element 
between the states of M = ±1/2. Thus the system of the type of Model I can not 
adiabatically change the sign of the magnetization. 
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